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Abstract 

We introduce a model for nonlinear viscoelastic solids where trav- 
eling shear waves with compact support are possible. We obtain an 
exact compact solution. We also derive a new Burger's type evolution 
equation associated with the introduced constitutive equation. 

PACS: 

94.05.Fg (Solitons and solitary waves), 

43. 35. Mr (Acoustics of viscoelastic materials), 

46.00.00 (Continuum mechanics of solids). 
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1 Introduction 



A compact wave is a non-linear solitary wave with a definite amplitude, which 
exists within the confines of a compact support; outside that support, it 
vanishes identically. Compact waves can be used to describe patterns with 
a compact support and sharp fronts. These are ubiquitous in Nature, but 
are hard to model mathematically. The driving force for the modeling of 
compact waves has been the 1993 discovery of compactons by Hyman and 
Rosenau [Ij. A compacton, in full analogy with the definition of a soliton, 
is a compact wave which preserves its shape and amplitude after a collision 
with another compact wave. 

Compact waves emerge essentially by mathematical degenerancy of the 
equations of motion, leading to a local loss of uniqueness, which enables 
the patching of two different solutions with a certain regularity. This point 
may be clarified and made rigorous in several ways, see Saccomandi [2] and 
Destrade et al. [3] who use the classical Weierstrass criterion to generate 
compact waves, and also Gaeta et al. [3], for example. 

From a purely mathematical point of view, there exist several equations 
that are factories of compact waves, such as for example the K(m, n) KdV 
equation [lj . On the other hand, there exist very few examples of physically- 
based equations which are capable of generating compact waves within the 
framework of a rigorous theory of material behavior. For solids, three ex- 
amples of such rigorous derivations have been presented so far: (i) in 1998 
Dusuel et al. [5j show that, in the continuum limit, the generalized $-four or 
double-well model with nonlinear coupling can exhibit compacton-like kink 
solutions when the nonlinear coupling between pendulums is dominant; using 
these ideas Saccomandi and Sgura [6] give an analysis of the full non-linear 
version of the classical models of DNA; (ii) Destrade and Saccomandi [7j use 
a non-linear theory of dispersion compatible with the axiomatic foundation 
of simple materials; (Hi) Goriely et al. [8] extend the theory of rods to non- 
linear material laws. In fluids the emergence o compact coherent structures, 
mainly via asymptotic method, has been a more investigated topic [9] . 

Here, we show that it is possible to generate compact waves in the frame- 
work of the nonlinear theory of viscoelasticity. To this end we consider a 
special viscoelastic theory of incompressible isotropic solids for which the 
viscous part is that of the Navier-Stokes theory, with a shear dependent 
viscosity — a common assumption in non-Newtonian fluid mechanics. 

In this framework, for a special class of constitutive equations, we provide 
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a rigorous existence result, an exact (albeit implicit) exact solution and an 
approximation for the compact kink. Moreover, we derive a new general- 
ized Burger's equation as asymptotic reduction of the full equation for the 
propagation of shear waves. 



2 Governing equations 

We call X the position of a particle in the solid in (£> r ), the reference con- 
figuration, and x the position of that particle at time t in (£>), the current 
configuration. A motion of the body is the one-to-one mapping x such that 
x = %(X,t). The deformation gradient F and the left Cauchy-Green tensor 
B associated with this motion are 

F =%£> B = FFT > (!) 

respectively, and the strain-rate tensor is defined as D = | (^FF^ 1 + F~ T F T ^j 

where the superposed dot denotes the material time derivative. An incom- 
pressible solid can undergo only isochoric motions, and this internal con- 
straint translates mathematically as: det F — 1, tr D — 0, at all times. 

We are interested in viscoelastic materials of differential type, with Cauchy 
stress tensor 

T = -pi + 2&B - + 20D. (2) 

Here, p is the indeterminate Lagrange multiplier introduced by the incom- 
pressibility constraint, /3i and are the elastic response parameters, and v 
is the shear viscosity coefficient. In all generality, (3% = h) (i = — 1, 1), 
where Ji, I2 are the first two principal invariants of the Cauchy-Green strain: 
Ii = tr(I?) and I 2 = tr(S~ 1 ). We assume that = (D • D), and moreover, 
that v > 0, i.e. the model is dissipative. 

The momentum equations, in the absence of body forces, take the form 
div T = pdv/dt, where p is the mass density and v = dx/dt is the velocity. 
Our aim is to investigate what happens in the shearing motion, specifically, 
x = X + f(Z,t), y = Y, z = Z, where the function f is as yet unknown. 
Straightforward computations give the components as 

/ 1 + K 2 K 

(By) = 10 

\ K 1 
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and 

( ° 

2 (Aj) = o 

where K = fz is the amount of shear, and the subscript denotes partial 
differentiation. Clearly we have now Ii = J 2 = 3 + K 2 and D • D = K t 2 /2. 

With ^ = 0, two of the three equations of motion are identically satis- 
fied. The remaining equation is pfu = dTi^/dZ. Therefore the determining 
equation for the amount of shear K becomes 

P K tt = \q(K 2 )K + u(K 2 )K t ] , (3) 

J ZZ 

where Q = 2(/3i + is the generalized shear modulus. The mathematical 
theory of quasilinear equations for viscoelasticity of strain-rate type can be 
found in ref. |10j . 

In order to rewrite eq. ^ in a dimensionless form, we need a characteristic 
frequency Q(> 0) so that it is possible to introduce the dimensionless time 
r = Qt; we also need a characteristic length L so that we can introduce the 
dimensionless length ( = Z/L. Usually the length L is determined by the 
geometry of the problem (e.g. the thickness of a slab wherein the wave is 
propagating). The characteristic frequency Q may be introduced in several 
ways: via the boundary conditions, by defining the ratio Q = Po/Oq, where 
/io = lim^-2_ s>0 Q is the infinitesimal shear modulus and z> = linift- t _>o or 
by a characteristic (finite) time t* at which, for example, localization of the 
solution occurs. Eventually, eq. ^ becomes 

SK rr = [QK + uK T ] cc , (4) 

where 5 = pQ 2 /fi L 2 , Q = Q/fio, and v = Q0/fi . 

As the final step, we specialize the equation of motion (via the constitutive 
relations) to the case of fourth-order elasticity. In particular, we take Q = 
1 + piK 2 , where p\ is a constant (/ii > for strain-stiffening solids and 
\i\ < for strain-softening solids), and assume the simplest form of shear 
viscosity dependence, namely, v = u + v x K 2 , where i>\ > is a constant. 
Note that pi = (/i + A/2 + D)/p, where p, A, and D are the second-, third-, 
and fourth-order constants of weakly non-linear elasticity [TTJ [TJ]. Under 
these assumptions, eq. Q reduces to 

5K TT =[K + fi.K 3 + u Q K T + Vl Kl] ?c . (5) 
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The existence and regularity of solutions for the Cauchy problem of such 
kind of equations has been considered by Friedman and Necas [TT]. Also, 
Pucci and Saccomandi p2] considered the quasistatic limit of eq. ^ and 
studied the mathematical and mechanical properties of the classical creep 
and recovery experiment. 

Eq. ^ is neither integrable nor linearizable, as may be checked via sym- 
metry arguments p3] so that we do not expect our solitary compact waves 
to preserve their shape and amplitude after colliding. In all likelihood, they 
are not compactons. 



In this section, we seek travelling wave solutions (TWS)s of eq. ([5]) in the 
form of kinks; i.e., continuous, bounded, monotonic waveforms that tend to 
constant, but unequal, limits at ±00. It is known that kinks may propagate 
in a viscoelastic medium; Jordan and Puri [T^l [T5] give an explicit char- 
acterization of such waves, and a detailed survey of the various qualitative 
properties of travelling waves solutions in viscoelasticity may be found in [16] . 

We begin our search for TWSs with the following observe: since eq. ^ is 
invariant under the transformation ( 1— > we need only consider, without 
loss of generality, right-traveling waves, i.e., solutions of the specific form 
K((,t) = <?(£), where g is a function of the single variable £ = ( — ct and 
the positive constant c denotes the wave speed. Substitution of the travelling 
wave ansatz into eq. ^ results in the following nonlinear ordinary differential 
equation (ODE): 



where primes denote differentiation with respect to £. Integrating this ODE 
twice, setting the first integration constant to zero, and then enforcing the 
usual (kink) asymptotic conditional yields 




3 Kinks 



(1 - 5c v + - w - "icwr = 0, 



(6) 





(1 - 5c 2 )g + ^(g 3 + g\g 2 + g\g x ) 



(7) 



Here, v = c~ 2 z/ /z/i and the wave speed is given by 





That is, g — > 31,2, as ^ — > T°°i where g± > g2 > are constants. 
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Having obtained the first order ODE satisfied by our TWS and deter- 
mined the wave speed, it is instructive to now consider the cases v = and 
v > separately. 

3.1 The case i/q = 

Focusing on this, the simplest case first, we further simplify the analysis by 
taking g 2 = and setting / = g/gi, in particular, c reduces to ci, where 



Cl = ^(l + ^gf)/5, (9) 

and eq. ([7j) becomes 

(ff = -af(l-f), fe[o,l), (10) 



where we have set a = c 1 3 /ii/z/ 1 for convenience. Now a standard stability 



analysis of eq. (10) reveals that / = 0, 1, the equilibrium solutions relevant to 
our investigation, are stable and unstable, respectively, for a > 0. Henceforth 
limiting our attention to only those dispersive solids that stiffen in shear (i.e., 



those for which fix > 0), we separate variables in eq. (10) and integrate. We 
are thus led to consider the quadrature 



/ 



<!/ -^ + /C, (11) 



[ /( l_ /2)] l/3 



where the integration constant /C will be chosen so that the kink is centered 
at f(0) = 1/2. Because of the zeros at / = and / = 1 in the denominator 



of the integrand, the left-hand side of eq. (11) is, in fact, a generalized inte- 
gral (see [3]). Nevertheless, it can be evaluated exactly in terms of special 
functions. Omitting the detail, it is readily established that the exact, albeit 
implicit, solution is given by 

2 iq(i, |; |; f)f^ - 2 iq(§, |; |; = -§^, (12) 

for £ e (£i,£o)> while outside this interval we have 



/(£) = {' ~; (13) 
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f 



tr= 0.25 




Figure 1: / vs. £ for cr = 0.25 and z/ = 0. Solid curve: eq. (12) and eq. (13). 



dashes): eq. (17)3. 



Broken curve (long dashes): eq. (17)i. Dots: eq. (17>2- Broken curve (short 



Here, 2 Fi denotes the Gauss hypergeometric series and the constants £i, £ 
are defined by 



-2^ + 3£i /4 & f ) 
6o-V3 



, £o = |<7 1/3 fii/ 4 (|, |), 



where 



B q (a,b) 



denotes the incomplete beta function. 

The shock layer thickness £ has the value 



lim /(£)- lim /(0 

5— >— 00 >+oo 



/|/(0)|=2(3a)- 1 /3. 



(14) 
(15) 



(16) 



If we expand the first term on the left-hand side of eq. (12) about / 



1,|,0 and then neglect the appropriate higher order terms, the resulting 
expressions can be solved for / in terms of £. Omitting the details, it is a 
relatively straightforward task, using the Heaviside unit step function, H(-), 
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to construct the respective approximations (see eq. (17)): 



/(e) 



#(£i-£) 



§(l-2£/*), 



1 - ^ [2nV3 - 



-1 2> 
3' 3/ 



3/2 



1^1 « I 



1 1. 4 

-3' 3' 3' 4 



4X4* 



1/3 



3/2 



(17) 

1 makes clear, the relatively simple approximate expressions given in 
eq. (17) are in very good agreement with the exact kink solution within, and 



As fig. 



even outside of, their theoretical ranges of validity. 

Remark 1: From eq. (16) it is evident that a shock wave, i.e., a propa- 
gating jump in /, forms as a — > oo since I — > in this limit (see fig. 2). In 
contrast, our compact kink solution does not exhibit acceleration waves, also 
known as "weak discontinuities" [TS1 §89], at £ = £ 10 because /'(£) £ C(J£); 
however, it should be noted that max|/'| — >■ oo as a — >■ oo (again, see fig. 2). 



3.2 The case z/ > 

Assuming now z/ > 0, we return to eq. Q and regard this ODE as a cubic 
polynomial in g'. Using Cardano's formula, the three roots of this cubic, 
each of which is a function of the single variable g, are readily determined. 
Fortunately, however, the cubic discriminant, which we denote here as D, is 
always positive. Hence, only one of the roots is real-valued. Denoting this 
particular root by Q(g), eq. becomes 

gf = g(g) = fe/2 + VD = (18) 

3^G/2 + 

where D = (^/3) 3 + (G/2) 2 and G = G(g) denotes the right-hand side of 
eq. ([7]). We observe here that Q{g) is strictly negative for g 6 (<? 2 , #1), by 
Descartes' rule of signs, and that g* = {—(gi + #2), #2, <7i} are the roots of 
both G(g*) = and £/(<?*) = 0, where —((71 + ^2) is an extraneous root in 
the present context. It should also be noted that the compact kink results 
discussed in the previous subsection are recovered with little difficulty by 
letting v — > (i.e., u — > 0). 
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-1 -0.5 0.5 1 



Figure 2(a-c): Plotted for a = 0.025,0.25,2.5, respectively, and 
Solid: — /' vs. £. Broken: / vs. £. 
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For general values of v, it is possible to show the existence of a kink 
solution, but it does not appear possible to determine its exact analytical 
representation. On the other hand, because Q = Gjv — G 3 /i) A + • • • , when 
expanded about G = 0, it is clear that as v — > 0, the points where G = are 
associated to a g' which is more and more vertical. This means that the kink 
compactifies as v — > 0. In other words, as the nonlinear viscoelastic part in 
the constitutive function becomes more important with respect to the linear 
part, the tails of the kink are of less importance. 

The process just described is clearly illustrated in fig. 3, where the param- 
eter A = vjg\ has been introduced for convenience. The sequence presented, 



which was generated from Eq. (18) by once again taking g 2 = and setting 
g — gxf, depicts the formation of the compact kink travelling wave profile 
(see Fig. 1) as v — > 0. 

Continuing under the assumption gi = 0, and with g = g\f, let us now 
expand Q = Q(f) about / = 0. On making the additional assumption A 



(j 2 / 3 and neglecting terms C(/ 4 ), eq. (18) is reduced to the following special 



case of Abel's equation: /' « —a\ 1 (f — f 3 ). Then, taking /(0) = 1/2, the 
"exact" solution of this ODE is easily found, using [TH eqs. (12)-(18)], to be 

1 + 3exp(2o-£/A) 



with a shock layer thickness of ^Abei = |A/o" that can never go to zero. Also, 



we observe that the graph of eq. ( 19 ) is very similar, qualitatively speaking, 
to the one shown here in fig. 3(a). 



Remark 2: From eq. (19) we find that, for sufficiently large values of A, 
and a zero limit as £ — > +oo, the viscoelastic model considered in the present 
Letter behaves very much like the cubically perturbed Kelvin- Voigt model; 
see, e.g., Ref. [Hj and those therein. 



4 A new Burger's equation 

Let us return to eq. ^ and introduce the new independent variables 9 = 
r — a C an d s = e£, where a = 5~ 1 ^ 2 and e is a small parameter. If we now set 
K{C^ t) = e 1 / 2 /«(s, 9) for some function k, then eq. ^ can be approximated 
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-4-2 2 4 



Figure 3(a-c): / vs. £ for A = 0.25,0.05,0.001, respectively, and a = 0.25. 
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via a generalized form of the modified Burger's equation (MBE)[^J namely, 

k s = p,i(K 3 ) e + VQK e e + i>i(Kg)e, (20) 

with the latter reducing to the former in the limit Ui — > 0. Here, we have 
set fii = yUia /2, z/ = z/ a /2, and v\ = via /2; and we have also assumed 
that u = 0(e). Additionally, we note that when u = 0, the structure of the 



travelling wave solutions of eqs. ( |5[) a nd (20) is exactly the same. We should 
also point out that, although eq. (]20[) is a second order differential equation, 
the compact kink is not an acceleration wave, because the jump in the second 
derivative is not finite. 



Moreover, because eq. (20) is an evolution equation, which is simpler than 
the wave equation given in eq. ([5]), other reductions to ordinary differential 
equations are easily found. For example, when v = 0, eq. (20) admits 



solutions in the separable form k(0, s) = ip(9)(j)(s), say. If uq = and /ii = 0, 



then eqn. ( 20 ) is a degenerate diffusion equation and the separable solutions 



are quite simple, with an interesting structure. Indeed, we have 

K(M= (21) 
^27(73 - S) 

Here, 7 is the separation constant and 71, 72, 73 are integration constants 
such that 247^1 —7 = 0, and therefore 7,73 > 0. Clearly, these solutions 
have a sharp front s = 72; and they blow-up in space for s = 73. 

It is interesting to consider the case of harmonic excitation and the cor- 
responding third-harmonic order generation. Given the initial condition 
k(0, 6) = k cos(u9), where k is a positive constant, we assume the solution 
is given by the sum of the fundamental and third-harmonic components; i.e., 
k = K\ + K3, where 

Ki(s, 6) = Re{ki(s) exp(iw6 1 )} = |/ti(s) exp(iu9) + c.c, 
ft3(s, 6) = Re{R,3(s) exp(3iu9)} = ^k 3 (s) exp(3iw6 ) ) + c.c. 

Here, "c.c." denotes the complex conjugate of the proceeding term and we 
require that \k 3 \ By successive approximation, we have 

^1 + ak x = 0, (22) 
as 



2 So named by Lee-Bapty and Crighton in 1987; see Refs. [El [14] and those therein 
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where a = u> u , and 

^ + 9cm 3 = -f^(i/ii-^i), (23) 

where the latter is solved subject to ^3(0, 6) = 0. 

Solving these ODEs in sequence using one of the many standard methods, 
we find, in turn, that 

Ki = Koe~ as cos(u9) , 
8a 

x \p lU sin(3u;fl) + z/^ 4 cos(3wfl)] . (24) 

These solutions reveal that the experimental measurement of the third har- 
monic in soft solids with shear-dependent viscosity gives direct access to the 
non-linear shear wave elastic parameter /2i and non-linear dissipation pa- 
rameter V\. The former is measured in the low-frequency regime, whilst the 
latter becomes dominant in the high-frequency regime. 

Remark 3: If we let v\ — > 0, and make the associations Ji\ h-> — |c~ 3 /3 
and vq 1 — y 5, then K\ and K3 respectively reduce to v\ and V3, which correspond 
to the MBE, given in pH eq. (42)]. 

Remark 4: In attempting to satisfy \ks\ <C [aci|, it is helpful to know 
that 

maxl/fal < fviuj 4 + \ui\u\ , (25) 

where 

Kq|z/icj 4 cos(3ws) + fiicj sin(3o;s) | 



max K3 = K3 

6»>0 



(26) 



and where 9* = (6a) _1 ln(3). 

Remark 5: The inequality max|K3| < sup|«i|(= Kq) is satisfied for Vw G 

(0, uj + ), where u + is the only positive root of uiu 4 + \p.\\uJ — 12kq 2 a\/2> = 0. 



5 Concluding remarks 

We have provided a clear constitutive framework that originates fields equa- 
tions admitting compact kinks. The constitutive framework is relevant be- 
cause it is a natural model for nonlinear viscoelasticity. The mechanism 
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of compactification is clearly related to the nonlinear viscosity term. We 
have also derived a new Burger's equation which maintains the compactifi- 
cation features of the full model. The usual approximation used in nonlinear 
acoustics to study the generation of the third harmonic indicates that the 
important difference between our model and the classical one is mainly at 
the higher harmonics. 

We established an important and original example of a how compactifi- 
cation may arise in the modeling of real-world phenomena. We point out 
that the mechanism generating compact waves here is completely different 
from those presented in P El E], where the compactification is made possi- 
ble by an interplay between nonlinearity and dispersion, and in j3], where 
compactification is made possible by considering non-smooth potentials. 
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